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Abstract. We consider the viscous n-dimensional Camassa-Holm equations, 
with n = 2, 3,4 in the whole space. We establish existence and regularity 
of the solutions and study the large time behavior of the solutions in several 
Sobolev spaces. We first show that if the data is only in then the solution 
decays without a rate and that this is the best that can be expected for data in 
. For solutions with data in n we obtain decay at an algebraic rate 
which is optimal in the sense that it coincides with the rate of the underlying 
linear part. 

Quelques questions de decroissance et existence pour les equa- 
tions visqueuses de Camassa-Holm. 

Resume : On considere les equations visqueuses de Camassa-Holm dans R", 
n = 2, 3, 4. Nous etablissons I'existence et regularite des solutuions. Nous 
etudions le comportament asymptotique des solutions dans plusieurs espaces 
de Sobolev quand le temps tend vers I'infini. On montre que si la donnee est 
seulement dans la solution decroit vers zero, mais la decroissance ne peux 
etre uniforme. Pour les solutions avec de donnee dans n f^™ on obtient 
une decroissance algebrique avec une vitesse qui est optimale dans le sens que 
c'est la meme que pour les solutions correspondant a I'equation lincaire. 



1. Introduction 
The Viscous Camassa-Holm equations (VCHE) are commonly written 
(1.1) vt + u-Vv + V -Vu^ + Vtt = v/\v 

u ~ a^Au — V 
V • u = 

Here we adopt the notation [v ■ Wu^)i = ^jVjdiUj. These equations rose from 
work on shallow water equations |3j , which led to [TT] , [T^] , where the equations are 
derived by considering variational principles and Lagrangian averaging. In light 
of this derivation the equations are sometimes called the Lagrangian Averaged 
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Navier-Stokes equations. In [9j, the equations were derived as a "filtered" Navier- 
Stokes equation, which obeys a modified Kelvin circulation theorem along filtered 
velocities. In this setting they are sometimes referred to as the Navier-Stokes- 
a equations, where a is the parameter in the filter. Solutions to the VCHE are 
closely related to solutions of the famous Navier-Stokes equation (NSE), but the 
filter allows bounds that are currently unobtainable for the NSE, making them in 
some ways better suited for computational turbulence study, see jl2j . 

In [S], [TD] these equations were studied in relation to turbulence theory, this 
treatment includes existence and uniqueness theorems on the torus in three dimen- 
sions. The two dimensional case was considered on the torus and the sphere in ^14] . 
Global existence and uniqueness in three dimensions was proved on bounded do- 
mains with zero (non-slip) boundary conditions in |^16J. These equations have also 
been studied in terms of large eddy simulation and turbulent pipe flow in [4], [5], [6], 
and [5] . In this paper we extend the current existence theorems and study the large 
time behavior of solutions. 

This paper is organized as follows. Section two consists of notation and conven- 
tions used throughout. Section three contains preliminary discussion of the VCHE 
and several useful lemmas. In section four we state existence and uniqueness re- 
sults for the VCHE, proofs of these statements are contained in the appendix. In the 
next two sections we continue the decay program of M. E. Schonbek, [5D], [3T], [25] . 
[25], [27]. The main result of chapter five considers solutions of the VCHE in the 
whole space and we prove that the energy of a solution corresponding to data only 
in i^(R") decays to zero following the arguments in 18j. We then demonstrate, 
by constructing counter examples, that no uniform rate of decay can exists which 
depends only on the initial energy. In chapter six we consider decay for solutions 
with initial initial data in n L^. We show, using the Fourier Splitting Method, 
that the energy of a solution decays at the rate expected from the linear part, this 
is the same rate of decay as solutions to the NSE. For solutions with initial data 
in iJ™ n we calculate the decay of derivatives using again the Fourier Splitting 
Method with an inductive argument. In section seven we examine how solutions of 
the VCHE approach solutions of the NSE strongly on intervals of regularity for the 
NSE. 



2. Notation 

In this paper, denotes the standard Lebesgue space with norm \\(p\\p = 
{J\4>\^Y^^- We use < u,v >= J uv to denote the standard inner product on 
the Hilbert space L^. Compactly supported solenoidal vector fields (subsets of 
Y, — {(p £ C^(ri)|V(/) — 0}) will be needed to describe incompressible solutions 
with zero boundary conditions. L^. will denote the completion of S in the norm 
II . lip. W"^'P will be used to denote the standard Sobolev spaces with the convention 
that i?" = (and = H°). The completion of S under the ff™ norm wiU 

be denoted by and {H^)' will be the dual space. To denote the Fourier Trans- 
form of a function (f> we will use either (f> or J-{4>), with or ^^^^(0) the inverse 
transform. Throughout we will use C to denote an arbitrary constant which may 
change line to line, to emphasis the dependence of a constant on a number, say 
we will write C(i^). 
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3. Preliminaries 

The Kelvin-filtered Navier-Stokes equations (KFNSE) are given by the formula 

— — h u ■ Vv + V ■ \/u + Vtt ~ vl\v 
at 

V -v = V ■u = 
v^Ou 

In the above, u = g * v represents a spatially filtered fluid velocity and O is the 
inverse of this convolution. The term u ■ Vw is similar to "mollifying" the Navier- 
Stokes equations, originally done by Leray, [15], to approximate solutions. The 
term v ■ \/u^ = ^Vj\7uj allows the solution to obey a modification of the Kelvin 
circulation theorem where circulation is conserved around a loop moving with the 
filtered velocity u. In two and three dimensions, using the identity 

(3.1) u-Vv + Y^ VjWuj = —u X (V X t;) + V(f • u) 

and including the term V(w • u) in the pressure, the KFNSE can be written as 

dv 

— + Vtt = u X (V X w) -f vAv 

V • u = V • u = 

v^Ou 

The following lemma will show that the bilinear term in the Kelvin-filtered Navier- 
Stokes equations behaves similar to the bilinear term in the Navier-Stokes equations. 

Lemma 3.1. Let u and v be smooth divergence free functions with compact support, 
then 

< u ■ \/v, u> + < V ■ Vu'^, u >= 
< u X {V X v),u >= 

Proof. The second equality is a consequence of the first, the identity (|3.1[) . and the 
fact that u is divergence free. To see the first inequality we just need to rearrange 
the terms and then integrate by parts 



/ VjdiUjUidx — —'Y^ I 



UidiVjUj dx 



□ 



Using this lemma, we can formally multiply the KFNSE by u to find 

. ^ d 

(3.2) < —V, u> +iy < Vw, Vu >= 

ot 

By choosing O to be the Helmholtz operator O = 1 — A we recover the Viscous 
Camassa-Holm equations 

Vt + u ■ Vt; + V ■ Vu"^ -|- Vtt = i^Av 
u — a^Au — V 
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In the case of the VCHE, (l3?2l) becomes 



(3.3) ^^(< u,u> < Vu, Vu >) + iy{< Vu, Vu > < Au, Au >) ^ 
This relation gives a priori estimates on u: 

\\ui;t)\\l + a^\\Vu{;t)\\l + 2iy f \\Vu{;t)\\ldt + 2iya^ f \\V^ui;t)\\ldt 

Jo Jo 

(3.4) < lluoll^ + a^llVuoll^ 

4. Existence of Solutions for the VCHE 

Existence and uniqueness of solutions for the VCHE on periodic domains in 
three dimensions was proved first in |10] using the Galerkin method. The most 
general existence and uniqueness theorems in three dimensions are provided in |16j 
which relies on a fixed point argument. The theorems in [T6] assume the initial 
data uq G Hq n with s G [3,5) and u = Au = on the boundary, where A is 
the Stokes operator. Here we state extended results which cover the whole space in 
dimensions 2 < n < 4, proofs are included in the appendix. As an intermediate step, 
we provide a new existence proof on bounded domains in dimensions 2 < n < 4, 
using the Galerkin Method, with initial data vq G L^, and u = w = on the 
boundary. Our bounded result in three dimensions is slightly stronger then [16j , by 
assuming vq G we have only implied u G H^. 

Definition 4.1. Let fl C M" be any open bounded subset or H. — R", n = 2, 3, 4. 
A weak solution to the VCHE (|1.1[) , with zero (no-slip) boundary conditions in the 
case of n bounded, is a pair of functions, u, v, such that 

vGL^{[o,nLl{n))nL'{[o,T]-H',{n)) 

u G L-([0, T]-Hlm n L2([0, T]-Hlm 

as well as v{x,0) = vq, and for any (j) G L^{[0,T];Hl{Vl)) with (j){T) = the 
following equalities are satisfied: 
T i-T 

< V, dt(f> > ds + < u - Vv, 4>> ds 
Jo 

i-T 



/ < 4>- Vu, V > ds + iy < Vw, > ds =< wq, ^^'(O) > 
Jo Jo 



and for a.e. t G [0, T], 



Theorem 4.2. Let Q C M" be an open bounded set with smooth boundary or 
— M", n — 2,3,4. Given initial data vq G H^^{Vl), AI > 0, there exists a unique 
weak solution to the VCHE il.l]) in the sense of Definition \4-l\ This solutions 
satisfies the estimate i3.4\ l as well 



as 



(4.1) |!9fV"^;(^)||^ + ^ [ \\dfV^+M^)\\l ds < C(||«o||hm) 

for all m + 2p<M. 
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Proof. Existence is given by Theorems 18.41 and 18. 61 in the appendix. The regularity 
statement is Theorem 18.81 and uniqueness is Theorem 18.91 The proofs foUow from 
the construction of approximate solutions using the Galerkin method on bounded 
domains. A priori bounds are obtained through energy methods. Using a com- 
pactness lemma we are able to find a strongly convergent subsequence which allows 
the limit of the approximate solutions to pass through the non-linearity. To extend 
to unbounded domains we solve the problem in balls of radius {Ri} (a sequence 
tending to infinity), and then invoke a diagonal argument. Regularity is established 
through an inductive argument relying on energy methods. □ 

Next, we will state a Corollary that describes the action of the filter and will be 
used many times in the following two sections. 

Corollary 4.3. 

\\dfS/"^u\\l + 2a2||arV™+iu||2 + a^\\dfV"^+Ml = Wdf^'^vWl 

\\dfn\\l + \\df\7^"+'u\\l < c\\d^v"^v\\l 

\\dfV"^uml + V f \\dfV"^+\is)\\lds < C{\\vo\\h^.) 
Jo 

for all m + 2k < M , where C is a constant which depends only on a,n, to, and k 
(in the last bound the constant depends also on ||?;o II 

Proof. This is an application of the Gagliardo-Nirenberg-Sobolev inequality to the 
bounds in the previous theorem. Differentiating the filter relation shows 

Squaring this relation then integrating by parts gives 

\\df\/"^u\\l + 2a'\\dfV^+'u\\l + a^dfV-^+Ml - l|5fV"Hl2 

This is the first bound in the corollary. Applying the Gagliardo-Nirenberg-Sobolev 
inequality to ||u||n and using the previous equality shows 

MV^uWl < C\\d!V^v\\l 

\\d!v"^+'u\\l < cwdfng 

This is the second set of bounds. Combining this with the regularity bounds in the 
theorem give the final set of bounds. □ 

5. Large Time Behavior of the VCHE: Non-Uniform Decay 

In bounded domains it is easy to see that the energy of a solution decays expo- 
nentially using the Poincare inequality 

\\u\\i<cm\^u\\i 

Indeed, start with the energy estimate p.3|) and apply the Poincare inequality to 
find 

i^(< u,u> +0^ < Vm, Vii >) + C{n)i^{< u,u> +c? < Vu,\7u >) < 
This differential inequality implies 

<u,u>W< Vm, >< C(||wo||2)e~^^"'''^* 
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The situation in unbounded domains is more delicate. If the initial data is assumed 
only in then the solution decays to zero but we are unable to determine the 
rate without more information, the precise statements of this idea are contained in 
Theorems [5?2] and 15.41 

We will follow [18j to show that the solutions in the whole space, constructed in 
Theorem l4.2l approach zero as time becomes large. The idea is to split the solution 
into low and high frequency parts using a cut-off function and generalized energy 
inequalities to show that both the high and low frequency terms approach zero. 
The idea of splitting into low and high frequency was first used in [17j . 



Lemma 5.1. Solutions of the VCHE constructed in Theorem \4-2\ with fl = R" 
satisfy the following generalized energy inequalities. Let E £ C^([0,oo)) and G 
CH[0,oo);Ci nL2(M"))^ then 



(5.1) E{t)\m) * v{t)\\l - E{s)U(s) * v{s)\\l + J E'iT)U{T) * v(T)\\ldr 

+ 2 [ E[t) < iP'{t) * v{t), iP{t) * v{t) > dr 

E{T)\\V^P{T)*v(T)\\ldT 

E{t) < u ■ Vu, iP{t) * iP{t) * v{t) > dr 
E{t) < V ■ Vu^ , '^{t) * 'iIj{t) * v{t) > dr 
For E £C^{[Q,oo)) and e C^{0,oo; L°°{R'')) we have 

(5.2) EmmmWl = E{s)ms)v{s)\\l + J' E'irmrMrmdT 

+ 2 f E{t) < tp'{T)v{T),tp{T)v{T) > dr 

J S 

-2 / E{T)Ui^{T)viT)\\ldT 

J S 

-2 [ E{t) < T{u-'^v),ij^{T)v{T) > dr 
E{t) < T{v ■ Vu^),iP^{T)v{T) > dr 



Proof. The proof of the first inequality is accomplished by multiplying the VCHE 
by E(t)^ * t/j * z; then integrating by parts and in time. The second inequality is 
obtained by first taking the Fourier Transform of the VCHE, then multiplying by 
ip'^v and integrating. □ 

Theorem 5.2. Let v be the solution of the VCHE il.l]) constructed in Theorem 
with n = M" and vq e L^(M"), then 



(5.3) lim||t;(t)|l2 = 
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Proof. We work in frequency space . We split the energy into low and high frequency 
parts 

(5.4) ||w||2 < 11(^^)112 + 11(1- 0)«||2 

with (j) — e^l''! will be chosen below. To estimate the low frequency part of the 
energy, begin with the generalized energy estimate (jS.ip . Temporarily fix t then 
choose E = 1 (the constant function) and 



g-|5|^(t+l-r) 



Note that tp and are rapidly decreasing functions for t < < + 1. The relation 

■0' = l^p-i/) shows the third and fourth terms in (|5.ip add to zero. Note (f> = e~l'>l = 
ip{t) and apply the Plancherel Theorem to see 



||<^*(t)||^ < ||el«l'(*-^V«(s)|| 



(5.5) 



J s 



With Holder inequality, Young's inequality, and the Gagliardo-Nirenberg-Sobolev 
inequality we bound 

\<4>^*u-Vv, e2^(*~^)f (t) >\<U^*u- V'y||2||e2^(*^^)i;(T)||2 

<C\m\_^\\u\\_^\\VvUvh 

< C(0)||w||2||Vu||2||Vi;||2 

Similarly, 

I < * w . Vu^, e^^'-'-^^vir) > I < 110' * ^; • Vu^||2||e2^(*""'i^(r)||2 

<Cm\_^\\v\\_^\\VuUvh 

<Ci^)\\v\\2\\Vu\\2\\Vv\\2 

Using the triangle inequality. Holder's inequality, and (|4.ip in (|5.5p yields 

110^(011^ < ||el«l'(*-^)#(5)||^ 

1/2 / \ 1/2 



2C(0)||«o||2 ^ \\yu\\ldT^ ^ llVi^ll^dr 



As the first term on the RHS tends to zero, applying the limit t oo yields 

/ poo \ 1/2 / „oo 1/2 

limsup|10«(O||^<2C(0)|lz;o||2 / llVuH^dr / \\Vv\\ldT 



The bounds p.4p and (|4.ip show IjVuljj and ||Vw||2 are integrable on the positive 
real line, letting s —^ oo leaves 

(5.6) limsup||0i)(t)||^ 

t — ^oo 

To estimate the high frequency start with the generalized energy inequality (j5.2p 
and chose i/i = 1 - e^l^l' = 1-0. Let Bait) = : |^| < G{t)} where G{t) wiU be 
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selected later and use < u ■ Vw, v >= to replace -0^ by 1 — V'^ in the 5th term on 
the RHS of (lOI). 



(5.7) E{t)\\Mml<m\\Ms)\\l+ t E'ir) f \Mr)\'d^dT 
+ [\e'{t)-2E{t)G\t)) [ \ijv{r)?didT 

Js Jsgir) 

+ 2 E{t)\ < Vw + w Vu^),(1-V;^(t))w(t) > | dr 

+ 2 f E{t)\ < T{v ■Wu'^),v{t) > \dT 

J s 

We remark both {1 — ip'^) and (j) = J-~'^{l — tp'^) are rapidly decreasing functions. Us- 
ing again Holder's inequality and the Plancherel theorem, then Young's inequality 
and the Gagliardo-Nirenberg-Sobolev inequality allows 

\ < T{u-Vv + v ■ Vu^), (1 - ljj^{T))v{T) > I 

= !<(!- V?(t))J^(w -Vv + v Vu^), w(r) > | 

< WJ'-^l - ^'(r)) *{u-\/v + v \/u^)\\2\\v\\2 

<c\\i~ ^^JM^jVvh + ||«||^||Vu||2)IK'||2 

<C{mv\\2\\yuh\\Vv\\2 

Similarly use Holder's inequality with the Plancherel theorem, then the Gagliardo- 
Nirenberg-Sobolev inequality, and Corollary 14. 31 to bound 

I < • Vu^),w(r) > I < ||w • Vii^||2||w||2 

< C||t.||_^||V?/||„||«||2 

<C\\v\\2\\Vv\\l 

Choosing E{t) = (1 + and G^it) = /3/2(l + 1) in ^1}, so that E' - 2EG^ = 0, 
and taking /? > sufhciently large, leaves 

/3(1 + r)'3-i 



||(1-0MOII2< 77^11(1 -'^)*(^)ll2 



\il -q^)v{T)\'d^dT 

Bcir) 



+ C||^;o||2^* |^i^||V«||2(||VH|2 + W^u^) dr 

For ^ e BG(i) and t sufficiently large, ■(/; = |1 - 0| < Therefore |1 - 0p < 

/3^/4(l -|- 1)^ and the second term on the RHS can be bounded as 



Js 4 Ja(t) Js 



/33(1 + t)- 



<^\\Mlj\^+r)-^dT 

<?ii«oii^(i + sr' 
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Letting t ^ oo shows 

(5.8) lunsupUl-<l,)v{t)\\l<^\\vo\\l{l + s)-^ 



+ C\\voUj \\Vv\\ldT + J WVuWldr) 

The bounds (|3.4p and (|4.ip again show |1 Vw||2 and || Vu||2 are integrable on the real 
hne. Letting s — > cx) proves 

hm sup 11(1 - (l))v{t)\\l = 

t — ^oo 

Combining this with (|5.6p and the Planchcrel theorem completes this proof. □ 



Corollary 5.3. Let v be the solution of the VCHE (jl.ip constructed in Theorem l4.2l 
with n = M" corresponding to vq e H^{R"). Then 

1 /•* 

lim - / ||t;(T)||2dT = 

t^oo t Jo 

Proof. Given an e > we can choose a large s such that ||w|l2 < e for r > s, this 
follows directly from the previous theorem. Then 

1 \\v{r)h dr=\f^ \HT)hdT + \fj \Hr)hdT 

(5.9) < \ r \\vir)\\2dT + e-~' 



Note that e was chosen arbitrarily and let i — > cx) to finish the proof. □ 



We have shown that the energy of a solution to the VCHE will tend to zero as 
time becomes large, now we will provide a counter example to show that there is 
no uniform rate of decay based only on the initial energy of the system. This is 
analogous to a result proved in [24j . The idea is to take a family of initial data with 
a parameter e that have constant norm, but norms of higher derivatives of the 
initial data can be taken arbitrarily small by picking e sufficiently small. It is then 
possible to bound the higher derivative norms of the solution arbitrarily small by 
taking e small.. Combining this with the energy relation (jS.Sp allows us to place a 
lower bound on the energy of the solution which depends on e. By choosing e small 
we can guarantee that a solution will remain away from zero for any finite amount 
of time. 

Theorem 5.4. Letv be the solution of the VCHE constructed in Theorem \4.S\ 

with n^M." and vq S (R"). There exists no function G{t, ^) : R+ x M+ M+ 
with the following two properties: 

\\v\\2<Git,\\vo\\2) 
(5.10) lim G(t, /?) = V/3 

t — *oo 

Proof. Fix uo{x) to be any smooth function of compact support and write Uq{x) = 
e^^^uoiex). Let Vq = Wq— a^Aitg and v'^ the solution of the VCHE given by Theorem 
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corresponding to the initial data vo- Note ||uq||2 — \\uo\\2 and ||V™w5||2 
e™||Vuo||2 for all e > 0. Also, 

(5.11) IKIl2 = holl2+«'l|V<||^ + a^||A<||2 

^\\uo\\j + a\^Vuo\\l + a^e^\\Auo\\l 

and 

(5.12) = llVu^ll^ + a^Aut,\\l + a^llVA^^II^ 

= e^WVuoWl + a^e^ll AmoII^ + aV||VAuo|l^ 



From the two previous inequalities and Corollary 14.31 we obtain a constant C = 
C(||uo||/f3), such that for all e > 

(5.13) \\vT2<C 

\\VvT2<Ce' 

Multiply the VCHE p.ip by Aw% then integrating by parts yields 

^^l|Vw1l2 + i^W^^v'Wl =< ■ Av'> + < Av' ■ Vu\v' > 

Use the relation < u% Vu^u^ >= 0, the Holder inequality, Sobolev inequality, and 
then the Cauchy Inequality to see 

\<u' ■ V^;^ Aw' > I |(-1) < (Vu<=) • V^;^ V^;<^ > | 
<'^\\Av% + C\\Vu^\\l\\Vv% 

Similarly, 

I < Av' -Wu^v" > I < C||Aw'||2||Vu'^||„|lwi^ 



-2 

<^\\Av^l+C\\WvT2\\^urn 

Applied to dsn]): 

(5.14) ^I'l^^'ll^ + i'l^^'ll^ ^ C||Vt.'||^||Vu'|l^ 



By ([331) and Corollary [43l 

\\yu^\\ldt<\\u^,\\l + \\vut,\\l<H\\l 

This bound, combined with (j5.13|) and ()5.14|) yields 

W'^vll < ||V«§||2e^ll"oll2 < Ce^e^^' 
Again, apply Corollary 14.31 

This together with the energy estimate p.3p implies 

^|(ll"1l^ + «^l|V«1|^)>-Ce^ 
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or, 

(5.15) \\u% + a^WVu^Wl > + a'\\Vut,\\j - CeH 

= \\uo\\l + e^a^\\Wuo\\l-CeH 
> \\ua\\l - Ceh 

From this we can deduce that there is no function G{t,(3,), continuous and ap- 
proaching zero in t for each fixed /3, such that ||u||2 < G(t, ||uo||2)- If there was such 
a function, then at some it would satisfy the bound G{tQ, Hwolb) < ||wo||2/2. By 
choosing e sufficiently small in (|5.15p . i.e. < ||Mo||2/4Cto, we have found initial 
data with a solution which cannot satisfy this estimate. □ 



6. Large Time Behavior of the VCHE: Algebraic Decay 

Although there is no uniform rate of decay for solutions with data exclusively in 
L^, we now show that there is a uniform rate of decay depending on the and 
norm of the initial data. Theorem 16.101 contains the most general decay result in 
this section. 

There is a relation between the shape of the Fourier Transform of the initial 
data near the origin and the decay rate of a solution to a parabolic equation with 
this data. By requiring the initial data to be absolutely integrable (in L^) we are 
in turn requiring the Fourier Transform of the initial data to be bounded. Using 
the Fourier Splitting Method it will also be shown that solutions in the whole 
space decay algebraically in as i ^ oo for initial data in n i/*^, M > 0. 
The decay obtained is the same as for the linear part (the heat equation). Note 
that the initial conditions can be weakened to require only that vq € X where 
X — {vo\vo{t) < C(l + t)^^} where vo{t) is the solution of the heat equation 
with initial data vq. The decay rate will depend on the relation between (3 and 
the number of dimensions. For similar results corresponding to the Navier-Stokes 
equations see [22], [29], and [30] . 

The Fourier Splitting Method was originally applied to parabolic conservation 
laws in [22] , and later applied to the NSE in [23] . In [21] the decay rate was made 
sharp in dimension n > 2 through a bootstrap method and logarithmic decay was 
shown for n = 2. In [31] the decay rate for n = 2 was made sharp through a 
bootstrap argument involving the Gronwall inequality. In this section we combine 
ideas from all of these papers in a slightly different way which allows us to prove the 
optimal energy decay rate in dimensions n > 2 without appealing to a bootstrap 
argument although we still use a bootstrap argument to obtain decay rates for 
higher derivatives. This same argument is also applicable to the NSE. 

The first goal of this section is to obtain a decay rate for the filtered velocity 
u, which is accomplished by applying the Fourier Splitting Method to the natural 
energy relation (|3.3p . This decay rate is then used with an inductive argument to 
obtain decay rates for the unfiltered velocity v and all of its derivatives. We start 
by finding estimates on HwHoo- 

Lemma 6.1. Let v be the solution of the VCHE constructed in Theorem \4-.S\ 

with fl = W\ corresponding to vq e Lin L^{W"-). Then, 



\T{v)\<C 



1/2 . X 1/2 

i+( / ii"(^)ii2) / w^vmi 
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where the constant depends only on the initial data, the dimension of space, and 
the constants in the VCHE (but not a). 

Proof. Use the identity 

(6.1) \j{uivi) = ujVvi + viVuj 

i i i 

and write the Fourier transform of the solution J-{v) as 

Jo 

where 

(6.2) *(e, t)=^-T{-K + Yl UiVi) -T{u-Vv~u- \/v^) 

i 

We would like first to bound "if, in that direction we have the following estimate 
which relies on the bound ||J^(0)||oo < ll^lli ^nd Young's inequality 

\T{u -Vv-u- Vt;^)| < C||m||2||Vu||2 

Also, taking the divergence of the VCHE (jl.ip shows 

A(7r + y^ UiVi) = div(uVw — uVv^) 

i 

Using the estimate immediately above and the Fourier transform leaves 
|^.F(7r + ^u,«,)| <C||u||2||Vz;||2 

i 

Now we can bound the integrand 

\^{i,t)\<C\\uU^vh 
Now take the supremum over ^ of (j6.2p and apply the Cauchy-Schwartz Inequality: 

\Hv)\<\Hvo)\+c( f \\u{s)\\lds\' (f\\^v{s)\\lds ' 



The bound |^(uo)|oo < Ikolli finishes the proof. □ 

Theorem 6.2. Let v be the solution of the VCHE constructed in Theorem 

\4-'^ with ft — M", corresponding to vq L'^ C] L^(M"). The solution satisfies the 
"energy" decay rate 

v-udx= ||w||2 + a^l|VM||2 < C(i + l)""/2 

where the constant depends only on the initial data, the dimension of space, and 
the constants in the VCHE (but not a). 

Proof. The previous lemma, with the bound (|4.ip . yields 



1+ / \\u{s)\\ids 



(6.3) < C 

Now we begin work with the energy estimate p.3p . Using the Plancherel Theo- 
rem we rewrite it as 
d 



dt 



il + a'\e)u^dC + 2i. / \^\^{l + a^\e)u'd^ = 
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Let B{p) be the ball of radius p where — f'{t)/{2vf{t)), and / is a positive, 
increasing function to be specified later. To simplify our equations we write E"^ = 



-/ E^di + 2.p^ 



BC{p) 



or 

(6.4) 



d 
dt 



E^ d^ + 2vp^ / E'^d^< 2vp^ 



E^d^ 



B(p) 



Recall the relation between u and v, that is v ^ u — a^Au which has Fourier 
Transform ii — v/{l + a^|fp). Combining this with (|6.3[) we see 



\E\\l, = 



< C 



\u{s)\\lds 



With this bound we can estimate the integral on the right hand side of (|6.4p . 



d_ 
dt 



E^ d^ + 2vp^ / E^d£^< Cp 



2+n 



1 



m(s)||2 ds 



We now have a differential inequality which can be solved using the integrating 
factor / to find 



n/2 



1 



|u(s)||2 ds 



Choose / = + go that /'// = (71/2 + 1)/(1 + t)and integrate in time from 

to r. 

(l + r)"/2+i/ E\^,r)d^< [ E^C,0)dC + C r{l+ f \\u{s)\\lds)dt 

Note < £'2 d^, then using the Tonelli theorem we can bound the integral 
on the RHS as 

(1+ / \\u{s)\\lds)dt< / (1+ / \\u{s)\\lds)dt 



r [ E'{^,s)d^ds) 

Jo JR" 



< ril + 

which leaves 

(l + r)"/2+i/ E^{tr)d^< f E'itO)d( + Cr{l+ I \\u\\lds) 



<C{l + r)+Cr f f E^{^,s)d^ds 
Jo Jr" 



This is of the form 

(t><C{l + r)+C{r) f (t>{s){l + s)-''^^+^ ds 
Jo 

with = (1 + r)"/2+i J^^ E^{£_, r) d^. The GronwaU inequahty now shows 
(l + r)"/2+i / E\^,r)d^<C{l + r)eMCr f {1 + s^^^-Us) 
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For n> 2 the integral r JJ (1 + s) "/^ ^ ds is bounded independent of r. Applying 
the Plancherel theorem one more time finishes the proof. □ 

Next we work out of order and establish the decay rate for the homogeneous 
norm of v using a similar argument as the previous theorem. 

Theorem 6.3. Let v be the solution of the VCHE constructed in Theorem 

\4-^ with fl — corresponding to vq g n L^(R"). The solution satisfies the 
decay rate 

\\yv\\i < + 

where the constant depends only on the initial data, the dimension of space, and 
the constants in the VCHE (v, a). 



Proof. Multiply the VCHE by Au, use the identity (|6.ip . after recalling that A-y is 
divergence free use the Holder inequality to obtain 

~\\yv\\l + u\\tiv\\l<C\\uU\Vv\\^^\\^vh 

After using the Sobolev inequality, Corollarv l4.3[ and the previous theorem, this 
becomes 



~\\Vv\\l + ^\\Av\\l<C{l+t)-/^\\Av\\ 



We will now restrict ourselves to t large enough so that C(l + t) ^ < v/2, this 
implies 

d_ 

It' 

The next step is to apply the Fourier Splitting method as in the previous theorem. 
Let B{p) be the ball of radius p where = /'/(^/) ^iid / is a positive increasing 
function to be specified later, using the Plancherel theorem: 



iVull^ + i/||Aw||^ < 



^Mi\\l + '^p'Ui\\l<'^p' f 

»t Jb{p) 



Lemma 16.11 with Theorem 16.21 imply 



\v\' 



< C 



With this bound the previous line becomes 



1+ / (l + s)-"/2ds 



(1 



1 + 



'^/^di 



WvWids 



Set / = (1 + and use it as an integrating factor 
dt 



((1 + tr^'+'iml) <C 1+ (^J\l + ds^ ^* II V^;|| 



Ids 



Again, as in the previous theorem, integrate in time from to r, then use the 
Tonelli theorem and the Plancherel theorem 

ft 



(l + r)"/2+2||Vv||2 < C{l + r) 



(l + s)-"/2ds ) dt] I \\\/v{s)\\lds 



The Gronwall inequality now shows 

(l + r)"/2+2||Vi,||2 < C(l + r)e^ 
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Where 

A = 

Note A is finite, hence 



\\yv{r)\\l < C(l + r)-"/2~i 



□ 



Corollary 6.4. Let v be the solution of the VCHE constructed in Theorem 
with f7 = M", corresponding to vq e Ll n ^^(R"). Then, 

\Hv)\ < C 
\Hn)\<C 

where the constant depends only on the initial data, the dimension of space, and 
the constants in the VCHE. 

Proof. Combine Lemma with Theorems 16.21 and 16.31 □ 

Corollary 6.5. Let v be the solution of the VCHE constructed in Theorem l4.2l 
with = M", corresponding to € n ^^(R"). Then 

\\v\\l < C(i + 

where the constant depends only on the initial data, the dimension of space, and 
the constants in the VCHE. 



Proof. In Theorem 16.21 we have shown that 

(6.5) \\u\\l + a^Wu\\l<C{t + l)--/^ 

Differentiating the Helmholtz equation, then squaring it and integrating shows 

\\Wu\\l + 2a^WMl + o^'\\'^Ml^\\'^v\\l 
Combine this with Theorem 16. 3|, 

With dm we see 

Ml = Ml + "ia^VuWl + a^\\Au\\l < C{t + 1)-"/^ 

□ 

We now turn our attention to a more general situation involving the Fourier 
Splitting Method. This next theorem will be used in the remaining decay proofs. 

Theorem 6.6. Let || V™?i;(0)||2 < oo. Given an energy inequality of the form 
and the hound 

m.t)\<c{i+tf 

which holds for < ■jjj^^pjy, we can deduce the asymptotic behavior 



V™w||^ < C 
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Proof. We proceed directly with Fourier Splitting. Apply Plancherel's Theorem 
and break up the integral on the LHS. 



B{p) 



Choose, for some large d, 

p' 



Then, using the assumption for the bound on w and performing the integration on 
the RHS we have 

^((1 + tfU^'wlll) <C[{1 + t)-™-l+d+2/3-«/2 ^ _^ ^y+d 

Integration in time and another application of the Planchcrel theorem finishes the 
proof. □ 

For the first application of the above theorem we will compute the decay rate 
for all spacial derivatives for solutions of the VCHE. 

Theorem 6.7. Let v be the solution of the VCHE constructed in Theorem 

\4-S\ with Q, — K", corresponding to vq G H L^iW""). These solutions satisfy the 
following decay rate for all m < K 



Proof. The cases m — 0,1 are Theorems 16.51 and 16.31 respectively. To prove the 
remaining cases, we first find an inequality in a form suitable for Theorem l6.6( then 
using inductive arguments establish decay. Having previously established regularity 
of solutions, we proceed formally. Let M < K then multiply the VCHE p.ip by 
A^^u and integrate by parts to find 

^||V*^t;||2 + iy\\V''+'v\\l < Im,o + Jm,o 



where 



, rn 

m=0 



M-m^ > 



Using the Holder inequality, the Sobolev inequality, CoroUarv 14.31 and the Cauchy 
inequality we find 



M 



Im,o = cJ2 I|V"u||„||V*^-'"«||^||V*^+1i;||2 
<C|M|^|lV*^+H-||2 + C||Vt,||2||V^'^«|l^ 



M 

+ CJ2 \\^"'-'v\\l\\V''+'-"^v\\l + ^\\V''+'v\\l 

m=2 
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We treat the other term in a similar way. 

M-l 

M-l 

Together, this leaves 

(6.6) iw^^^'^Wl + < c\\v\\l\\^''^'v\\l 

+ C\\Vv\\l\\V''v\\l 

M 

+ Cj2\\S7^^-'v\\l\\S7^+'--v\\l 

m=2 

The remaining part of this proof will proceed by induction where the base case is 
Theorems 16.31 and 16. 51 We assume (inductive assumption) that the decay 

llV™w|l^ < C(t + l)-™-"/2 

holds for all m < M and will show that it holds for m — M. With this inductive 
assumption (|6.6p becomes 



(6.7) j^\\^''v\\l + ^\\^''+'v\\l < C{1 + t)-/^V^'+\\\l 

+ C(l + i)-l-"/2||V*^w||^ + C(l + i)-A^-"/2 

Consider t large enough so that C(l + f)^"/^ < Subtracting the first term on 
the RHS, (l6Jl) becomes 



^||V*^t;||2 + ^\\W^I+^v\\l < C{1 + t)-l-"/2||y 



A/„.l|2 



2 



dt" 4' 

+ c(i + t)-*^-"/2 

The next step is to apply the bound ||V^^w||2 < C (Theorem 14.21) with Theorem 
16.61 to obtain the decay rate 

||V^-^t;||^ < C(l +t)-"/2 

Continuing with a bootstrap argument, placing this new bound into (j6.7p and again 
using Thcorcm l6.6l the optimal decay rate is obtained and the proof is complete. □ 

The next goal is to extend the decay results to time derivatives of the solution. To 
begin we will compute a frequency bound for the spacial derivatives of solutions to 
the VCHE. This next lemma will be used inductively with Theorem l6.6l to compute 
decay rates for the norm of all time derivatives. 

Lemma 6.8. Let P > 1 and v be the solution of the VCHE U.l\) constructed in 
Theorem\4^wtth n^W, corresponding to vq e n Li(M"). // 

\\dy"'v\\j < C(l + t)-2p-'n-"/2 

for all p < P and m — 0,1, then 

\d^viO\<c{i + t)-^ 
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for < -jfj^- Here the constant depends only on the initial data, the dimension 
of space, and the constants in the VCHE. 

Remark 6.9. Note that the conclusion for P = is true by CoroUarv l6.4l 
Proof. The chain rule 

applied to (|6.2p shows 

+ f\-mYe'^*-^^\^\''^{^,s)ds 
Jo 

We bound 5* (defined by 16. 2p similar to the proof of Lemma 16.11 but using the 
assumptions of this Lemma. 

df^{tt)^dfA + dfB + dfC 

3 

1=0 

< C(l +t)"P""/2-l/2 

\d?B\ = \dfY.Tiu,VvJ)\ 

j 

<j2c\\divh\\dr'yvh 

1=0 

< 0(1 + typ-'^/^-^/^ 

\dfC\ = \d^^H^)\ 

< \dfA\ + \dfB\ 

< c{i + t)-p-'''/^-^/^ 

The bound < C (Corollary I6.4p and |C| < ^ finish the proof. □ 

Theorem 6.10. Let v be the solution of the VCHE U.l\) constructed in Theorem 
\4-^ with fl — M", corresponding to vq G L^(M."). These solutions satisfy the 
following decay rate for all m + 2p < K 

\\dy"'v\\i < C{t + l)-2p-™-"/2 

where the constant depends only on the initial data, the dimension of space, and 
the constants in the VCHE. 
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Proof. This proof follows closely the proof of Thcorem l6.71 we first find an inequality 
in a form suitable for Theorem 16.61 then using inductive arguments we establish 
decay. Choose P and M such that M + 2P < K, then apply df to the VCHE 
p.ip . multiply by dflS^'^v and integrate by parts to see 

^^Wd^V^'vWl + v\\dfv''+\\\l < Im.p + Jalp 



where 



D—O ni—0 
P M-1 



p ) \m 



p=0 m=0 ^ ^ 

or, in the case M — 0, 

p 



p=0 

Use the Holder inequality, the Sobolev inequality, Corollarv l4.31 and the Cauchy 
inequality we find, for M > 0, 

p M 

p^O m— 

< cE \\oMi\\dr'^''^'v\\i +cj2\\df^-\mdr'^''v\\i 

p=0 p=0 
P M 

p=0 m=2 

Similarly for the second term if M > 0, 

Jalp < ^E E ll^f V^^+i«|l2|iafV™+iu|l„||ar-^V^^-™t;|l ^ 

p=0 m=0 

p=Q p=0 
P AI-1 

+cE E ii5rv"-^.ii^ii5r-^v^^+i-'%ii^+^ii9rv^^+i.ii^ 

p=0 m=2 

In the case M = the estimate is 

p 



p=0 
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We have shown in the case M > 
(6.8) 

+ cf2\\df^^\\l\\dr''^''v\\l 

p=0 
P M 

^P-PYjM+l-m 



p=0 m=2 



V 



and in the case M = 0, 



p 



(6.9) + "^wdrvvwi < cJ2\\d?v\\i\\dr'yv\\i 

We now begin the inductive part of our argument where the base case is Theorem 
16.71 Pick P < K/2 and assume (inductive assumption) the decay 

(6.10) l!5fV™z;||2 < C{t + i)-2p-™-«/2 

holds for all p < P and m such that 2p + m < K. The inductive claim is that 
the decay holds for p = P with m such that 2P + m < K. To prove the inductive 
claim it will be shown first that the decay rate holds for p = P and m = using 
(j6.9p . Then, using (|6.8p it will be shown that the decay rate holds for the remaining 
values of m using another inductive argument. 

To establish the decay for p = P and m — 0, apply the inductive assumption to 
(EH) to find 

j^WoMl + "^WdfwvWl < c{i + t)-"/^drwv\\l + c{i + 1)-'^-'- 

Take t large enough so that C(l + i)^"/^ < and move the first term on the 
RHS to the left side 

j^\\d^'v\\l + ^\\drS7v\\l<C{l+tr'^-'-- 

Now, an application of Theorem 16.61 with Lemma 16.81 establishes the decay (jG.lOp 
for p = P and m = 0. This is the base case for the next inductive argument. 
Assume (inductive assumption) the decay (|6.10p holds for m < M + 1 when p < P, 
and m < M when p = P, we will show that this implies the decay holds for m = M 
and p = P. Proving this inductive claim will finish the proof. Begin by applying 
the inductive assumption to ([67 



^ l|2 I '^IISiPvtM+I ||2 ^ /^fi I ,\-ri/2|| QPyyM+l ||2 

— ||dt V f||2+2ll'^t^ ^ v\\2<C{l+t) ' lldf V ^ V\\2 

+ C(l+i)-"/2-l||aPvA./„||2 
+ C(l + t)-2p-A/-n 

Take t large so that C(l + i)^"/^ < i^/4 and move the first term on the RHS to the 
LHS. Then apply Theorem 16.61 with Lemma [6751 to establish the decay rate 

(6.11) ||9fV'"^;|l^<C(i + l)-"/2 
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Another bootstrap argument gives the optimal decay and finishes the proof. □ 

7. Convergence of the VCHE to the NSE in the Whole Space 

To understand how solutions of the VCHE approach solutions of the NSE as the 
filter constant a approaches zero we must first understand how a solution u of the 
Helmholtz equation 

(7.1) u-a'^Au = v 

approaches v when taking a to zero. To begin we state a theorem concerning the 
Helmholtz equation in all of space, the theorem is standard elliptic theory and no 
proof is given. This theorem can be proved using elliptic estimates and interpolation 
or if one multiplies the Helmholtz equation by e"^/" and divides by it can be 
thought of as the heat equation and the bounds follow from estimates on the heat 
kernel. 

Theorem 7.1. Given v G LP(R"), p e (l,oo), there exists au £ W^^'P(R") that is 
a weak solution to the Helmholtz equation u — a^Au = v. Moreover, this function 
satisfies 

\\u\\p < \\v\\p 

,, ,, C(n, v.q) ,, ,, „ n A 1 , 

M, < ^Jii^hh f or ^ =-{---)< I 

Ifn{2/p— 1) < 1 then the solution is unique. 

Proof. Standard elliptic theory. □ 

A solution u of the Helmholtz equation corresponding to v will approach v weakly 
as the filter parameter tends to zero. Indeed, fix w g iP(R") and let {a^} be 
a sequence tending to zero. By the above theorem, for each there is a weak 
solution G W^'P{M.") of the Helmholtz equation such that 



/ Ua- ■ (f)dx + af I S/uai ■ V0 dx — V ■ 

JR" JR" JR" 



The functions Uq,. are bounded in ^^(R") independent of a^, so there exists a 
(possible) subsequence ai- with a weak limit in LP(R"). Also, for l/p+ \/q = 1 

[ Vu-V<j)dx <a^\\Vu\\p\\V(^\\g<Cin)ay^\\v\\p\\V<j)\\q 



which approaches zero as 0. This proves that v in iP(R"). We can do 

better then this if v is sufficiently differentiable. 

Theorem 7.2. Let v £ M^^'^(R") and let u he the corresponding solution to the 
Helmholtz equation |7.i[ j. Then 

\\u -v\\,< C{n,p, q)a'/'-^\\Wv\\p /or 7 - ^(1 - 1) < 1 

2 p q 2 

If a is a sequence tending to zero and Uq, are solutions the Helmholtz equation, then 
Ua V strongly in L'^(R") for 1/p — 1/q < 1/n. 
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Proof. If u and v satisfy the Helmholtz equation, then 
(7.2) \\u-v\\g <a^\\Au\\g 

The Hchxiholtz equation is Hnear, so the derivatives of the functions obey the rela- 
tion Vm — a^AVu — Vw. Applying Theorem 17. II to this PDE with the restriction 
on 7 allows the bound 

II A II ^ C{n,p,q) 
\\^^h< ^3/2+7 II^^IIp 

Together with (|7?2l) . 

\\u-v\\,<C{n,p,q)a'^'-'\\^v\\p 
The second statement is an immediate consequence of this. □ 

In [S] , [TU] , the authors show how the solutions of the VCHE approach a solution 
of the NSE weakly when the parameter in the filter tends to zero {a ^ 0). We will 
show how solutions to the VCHE approach solutions to the NSE strongly as a ^ 
when the solution to the NSE is known to be regular. The proof requires estimates 
on the solution of the VCHE which are independent of a, but in regions of time 
where the NSE is known to be regular by some functional analytic arguments, the 
passive bound on the filter make this assumption reasonable. 

For example, solutions of the Navier-Stokes equation obey the Prodi Inequality 

m 

jjVu\\l<C\\Vu\\r 

This can be used to prove existence of a strong solution in some time interval [0, T] 
or regular solutions for all time if the initial data is small. The Prodi inequality is 
proved through energy estimates, using the passive bound for the filter in Theorem 
17.11 and following the same energy arguments allows the same bound for solutions 
of the VCHE. This bound will be independent of a, so we can apply the following 
theorem to conclude that in some closed interval [0, T] the solution of the VCHE 
approaches a solution to the NSE strongly. 

Theorem 7.3. Let {a^} be a sequence of filter coefficients tending to zero and 
the solutions of the VCHE il.l}) constructed in Theorem \4-S\ with Q = M" 
corresponding to wq G _ff^(R"). Let w be the solution the NSE with initial conditions 
wq. In any time interval [0, T] where a solution to the NSE is known to be regular, 
if there exists a bound 

sup sup {\\Va,\\l + WVVa^Wl) < C 
Oi tG[0,T] 

which is independent of a, then Va approaches w strongly in L°°([0, T], L'(R")) as 
a ^ 0, where q = 

Proof. We begin with a mild form of the solutions to both problems. We are 
working in a time domain with known regularity so these are the unique solutions. 
If P is the Leray projector onto the divergence free subspace of L^ and $ is the 
heat kernel, then 

w{t) = $(t) *wa- [ <I"(t - s) *¥[w Vw] (s) ds 
Jo 

v{t) = $(i) * Wo - / <I>(t - s) * P u-Vv + y2 UjVvj (s) ds 
Jo '- 
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By adding and subtracting cross terms we see 

w{t) - v{t) = - I <^{t - s) [{w ~ u) ■ Vw + u ■ V(w - v)] 

+ P [uj\'{vj — Wj) + {uj — Wj)'S/wj\ (s) ds 

The first term in the integrand is bounded using Young's inequahty and the defi- 
nition of the projector 

- s) * P [(w - m) • Vw] (s)||, < - s)||p|| {w-u)- Vw\\2 

< \\<^>{t- s)\\p\\w-umVw\\i 

where l/q + I = 1/p + 1/2 and 1/2 = l/q + l/l. Using Theorem O with 7 = 
(1/2 - l/q)n/2 < 1/2 we obtain 

\\<f{t-s)*F[{w - u) ■ Ww] < \\<i>{t^s)UVw\\i (\\w - vWg + Ca'/^-'WVvh) 

The fourth term can be bounded in essentially the same way. We approach the 
second term in a slightly different way, by first passing the derivative to the heat 
kernel. These functions are smooth functions of the whole space so the projector 
will commute with the derivative. 

||<I>(t -s)*P[u- W{w - v)] \\g = II V<I>(i -s)*F[u-{w- v)] ||, 

Then by using Young's inequality and the definition of the projector 

||$(i -s)^F[u-V{w- v)] \\g < ||V$(t - s)||p||w|Mh - V\\g 

To bound the third term, start with the product rule and again pass the derivative 
the heat kernel 

mt ~ s) * P [u.Viv, ~ w,)] II, = ||V$(t -s)^P [u,iw, - V,)] II, 

+ \\<l>{t-s)*F[Vu,{w,-Vj)]\\g 
Then, using Young's inequality and again the definition of the projector 
\Mt-s)*¥[u,Viv,~w,)]\\, 

< (||V$(i ~ S)||p||7.,||, + ||$(i - ■s)\\p\\VUj\\l) \\Wj - v,\\. 

Putting all of these bounds together and estimating the heat kernel yields 

||«-w||g < + B / \\v-w\\g{s)ds 

Jo — sj 

r* 

A = C I WWvhds 



B = C sup (||Vw||( + ||u||; + ||Vw||i) 
se[o.T] 

Here, S — 1/2 + (1 — l/p)n/2 < 1 by the assumption I > n. Application of the 
Gronwall inequality finishes the proof. For example, a modified Gronwall inequality 
|26j now shows 

||f ~w%< Aai/2-''T(Br(l - 6)t^) 

oc 

See also [I]. Letting a — * we see that v ^ w strongly in □ 
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8. Appendix 

Here we construct a weak solution to the VCHE. Due to the close relation be- 
tween the VCHE and the Navicr-Stokes equation, our proof is similar to known 
existence proofs for the NSE. See, for example, [1], [2], [7], [13], [15], [28]. First, we 
construct solutions on any bounded Q with smooth boundary using the Galerkin 
method, this is where the Stokes operator is known to be compact thanks to the 
Poincare Inequality. Special care is taken to use inequalities which do not depend 
on the size of fl so we can use these solutions to prove existence of a weak solution 
in unbounded domains. The only step that requires Q bounded is in the compact 
inclusion used to obtain the strong convergence necessary to pass limits through 
the non-linear term. This problem is overcome by working in the compact support 
of the test functions. 

To begin we recall a standard but useful elliptic estimate. 

Remark 8.1. Let 2 < n < 4 and fl C K." be an open set with smooth boundary. If 
u e H^, and v G L'^ satisfy the Helmholtz equation 

u — a^Au ~ V 

on Q, then 

(8.1) IM|„<C||t;||2 

(8.2) ||Vm|1„ < C||i;||2 

(8.3) + 2a^\\Vu\\l + a'^H Au||^ = 

where the constants C depend only on a and n. 
The stationary Stokes equation 

Am -|- Vp = V 
u\dn = 

is known to have a solution u £ H^{n) for each v S L^i^) when is an open 
bounded set. Solving this PDE defines an operator L'^(^) i?^(il). Composing 
this with the compact inclusion H^{fl) L'^i^) gives a compact and self-adjoint 
operator L'^{i^) ij(r2), which we call the Stokes operator. 

Lemma 8.2. Let C M" be an open bounded set. There exists an orthonormal 
basis of L^{fl), {wj}!^!, where each ujj is an eigenfunction of the Stokes Operator 
on fi. The associated eigenvalues are all positive real numbers and the eigenvectors 
are smooth and approach zero on the boundary. Let H„i — span{LOi, w,„} and let 
Pm be the orthogonal projection P„i : L'^{n) —^ H„i- Given vq G C^{fl), for each 
m there is an approximate solution 

m 

Vrn = y^^gjm{t)u!j 

i=i 

and 

E9im (t) 

7 = 1 J 
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where gjm G C^([0, Tm]) for some time Tm- These approximate solutions satisfy 
the following relations: 

(8.4) < dtVm, LUi > + < Um ■ Vw^, LUi > ~ < UJi ■ WVm, U„i > = U < Av 

Vjn{0) = PmVo 

Proof. Owing to spectral theory the Stokes operator (self-adjoint, compact) has a 
countable number of positive eigenvalues , and associated smooth, divergence- free 
eigenfunctions oji which form a basis for L^(il). These functions satisfy the relation 

- Auji = XiLOi 

To determine the scalars gim we construct a system of m ODE's. 

— r- h I'Aigim 

at 

m 

I \ ^ 9jm9km I y-, y-, 

+ 2^ -, ax ' ^^i'^* > - < • yujj,ujk >) = 

3,k=l 

Local existence of solutions to ODE's give existence of solutions gim, which are 
defined for some time interval [0,Tm]. □ 

The bounds in the next lemma will prove that T,„ can be bounded independent 
of TO, and in fact Tm = oo for all m. 



Lemma 8.3. For 2 < n < 4, the approximate solutions constructed in Lemma \8.2\ 
have the following bounds, which do not depend on T, f2 or m. 

\\vm\\L'^(lo,T]-Ll{n)) + II V^;m||L2([o,T];L2(o)) < C'{n, a, v, Hi^olb) 
||9tW™||L2([o,T]:(ffi)'(o)) < C{n,a,v, Hi^olb) 

Proof. Similar to formal multiplication of the VCHE by u, multiply (|8.4p by 
jqi^sy-ffim, sum, then apply Lemma l3.1l to see 

\\u^\\l + a^\\'^Um\\l + 2v / llVw^ll^dt 



,2. 



.5) +2a'v ||Aw„||^d< = ||mo||2 + « llVuol 



2 



With the Poincare inequality and (|8.3I this becomes first bound in the theorem. 
Using WY\ we deduce 



(8.6) lkm||„+ / ||Vu™||^,dt < C(n,a,z/,||uo||2,||Vuo||2) 

Jo 

To bound the derivative start with (|8.4p . Any e can be written as a sum 
of the uji so each approximate solution satisfies 

< dtVm,4> > + <Um - VVjn,4> > - < (f> ■ VVjn,Ujn >= V < AVm,4> > 

After integration by parts and applying the Holder inequality with the Gagliardo- 
Nirenberg-Sobolev inequality we find 

||2||V0||2+C||Vf™||2||V,^||2 

As (j) was chosen arbitrarily we conclude 

\\dtv,n\\(Hiy < C(||u„,||„||Vt;™||2 + ||V?;™||2) 
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This, together with (j8.5|) and (|8.6p . proves the second bound in the theorem. □ 

Theorem 8.4. Let Q e R", 2 < n < 4 6e a bounded set with smooth boundary and 
vq € C^(ri). Then, there exists a weak solution to the VCHE il.l]) in the sense of 
Definition CTTp. 



Proof. Thanks to Lemmas 18.21 and 1 8 . 31 we only need to prove the convergence of the 
approximate solutions. Lemma 18.31 shows how the sequence Vm remains bounded, 
so using a possible subsequence and the Banach-Alaoglu theorem there exists a 
function 

v€L°°i[0,T];Ll{n))nL'i[0,T],Hlin)) 
dtveL^{[0,T];(Hlnn)) 

such that 

(8.7) Vm in L^{[0,T]:Ll{n)) weak* 

(8.8) Vm^vin L'^{[0,T];Hl{n)) weakly 

We will now show that u is a weak solution to the VCHE p.ip . 

By the construction of our approximate solutions and integration by parts, we 
know for any basis vector Uj g L'^i^) and any smooth scalar function of time 0j(t) 
such that (j)j{T) = 0, 

.T .T .T 

/ < Vm, ip'^j > ds + < u ■ Vw, (jjLUj > ds + < (j)ujj ■ Vu, V > ds 
Jo Jo Jo 

< Vv,V(j)ujj > ds u,„(0), (/)(0)wj > 

The convergence (|8.7|) and (j8.8|) implies 

ft pt 

< Vm,(t>'j^j > ds ^ < V,4''jLLJj > ds 

Jo 

/ < S/vm, > ds ^ < \7v, (j)j\7ujj > ds 

Jo Jo 

Also, 

(8.9) < Vm{0), </'i(0)wj >=< Pm{vo), </'j(0)tJj >-^< vn, (j)j{0)ujj > 

Passing through the non-linear terms will require strong convergence, so we use 
the fact that the bounds in Lemma [8731 implv (see [7 , Lemma 8.2) the existence of 
a possible subsequence Vm such that 

(8.10) Vm^v in L^{[0,T];L^{n)) strongly 
Theorem 17.11 give the existence of a function u which satisfies 

u — a^Au — V 

Similar tolO 

\\um ~u\\l + a^||V(u™ - u)\\l + a^\\V'^{um - u)\\l = \\vm - v\\l 

In particular, applying the Gagliardo-Nirenberg-Sobolev Inequality shows \\um — 
^ C||wm — ^112- This, with the strong convergence (|8.10p . shows how Um 
approaches u strongly. 
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We can now prove the convergence of the non-hnear terms 

/ < u„i ■ Vwm, 4>jUJj > ds ^ < u ■ Vu, > ds 
Jo Jo 

t nt 

< cf)jUOj ■ S/Vm,Um > ds ~^ < (j)jUJj ■ Vw, U > ds 



Indeed, adding and subtracting the cross terms, then using the Holder Inequahty, 
the Gaghardo-Nirenberg-Sobolev Inequahty 

\ < Um ■ VVm,(t>LjJj > — < U ■ Vv, <j>jUJj > \ < Ai + Bi 
Ai = \ < (u,„ - u) ■ VVm, 4>jUJj > I 

< ll^m - u|l2||Vw„||2||'/)jVcJj||2 

Due to the strong convergence (|8.10p the bound in Lemma 18.31 and the Holder 
inequahty, we see Aids 0. Similarly, 

Bi = \ <u - V(u™ - v), (f>jUjj > I 

= \ <U - V0jWj, {Vjn -v) > \ 

< \\u\\n\\(f>jVl^j\\^\\Vrn - vh 

< C||u||2||0jVWj||_2zi_ ||w„i -v\\2 

Again, owing to (|8.10p . Lemma [531 and the Holder inequality, Bids 0. 
Putting this together, 

i-T pT 

< Um ■ VWm, (f^j^j > ds ^ I < U ■ Vw, (pjiOj > ds 

Jo 

The remaining non-linear term is handled in a similar way 

1 < (f)jU!j ■ Wv„i,Um > — < (l)jLLjj ■ Vu, U > \ < A2 + B2 

A2 = \ < (jjjUJj ■ VVm, [Um - w) > I 

< ||(?!'jC^j||^||Vu,„||2||Mm - U\\n 

< C\\(l)jVuJj\\2\\'VVm\\2\\Vrn " v\\2 

B2 = \ < (t)jujj ■ V(u™ -v),u> \ 
= I < (f)jujj ■ Vu, v„i-v>\ 

< \\(f>jUJj\\^\\Vm - v\\2\\yu\\n 

< C||(/)jVWj||2||Wm - 'y||2||w||2 

Applying (|8.10p with Lemma 18.31 and the Holder inequality shows 

' < (pjUJj ■ VVm, Um > ds ^ < 4>'jOJj ■ Vw, U > ds 

Jo 
Since the ujj are dense in and is an arbitrary smooth function the proof is 
complete. □ 

Corollary 8.5. The conclusions of Theorem 18.41 hold with the relaxed hypothesis 

vo e Llin). 
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Proof. Note that all of the bounds attained in Lemma 18.31 and used in the proof 
of the previous theorem depend only on the norm of the initial data. Let 
Vq e C^(ri) be a sequence of functions approaching vq strongly in Hq such that 

Such a sequence can be constructed using standard moUifiers and cutoff functions. 
Considering each Vq as initial data, Theorem l8.4l and its corollary give the existence 
of a weak solution in the sense of Definition 14.11 Applying (|8.3p . we see that 
these weak solutions satisfy the bounds 

\\'"'\\L-^(lo,T];Ll{n)) + II Vt;lL2([o .r].i2 (Q)) < C{n, a, ly, Hwolb) 

l|5fylL2([o,T];(Hi)'(o)) < C{n, a, iy,\\vo\\2) 

and for each cj) G the relation 

(8.11) 

i-T rT 

< v\ dt4> > ds + < • Vw*, (/) > ds 
Jo 

T 

< 4> ■ Vu',i;' > ds + < Vw', V(/) > ds vq, 4> > 



Jo 

As before, using the Banach-Alaoglu Theorem and extracting a possible subse- 
quence implies that there exists a function 

veL^{[0,nLl{n))nL\[0,T];Hlm 

dtv&L'i[0,T];{Hlnn)) 

such that 

v' ^vin L°"{[Q,T];Ll{n)) weak* 

v' ^vin L^{[0,T];Hl{n)) weakly 

Passing the limits through (|8.1ip follows by the same steps as in the proof of the 
previous theorem. □ 

Theorem 8.6. Let vq € L^(M") Then, there exists a weak solution in the sense of 
Definition\4-.l\ with initial data vq in the whole space K", 2 < n < A. 



Proof. Let Ri be a sequence tending to infinity and XRi ^ smooth cutoff function 
which is equal to 1 inside the ball of radius Ri — e and zero on the boundary of the 
ball with radius Ri . The Corollary 18.51 now gives existence of a weak solution v^^ 
on the ball of radius Ri with initial conditions v^XRi ■ Extend v^^ to all of M" by 
setting it equal to zero outside the ball of radius Ri. All of the bounds in Lemma 
8.31 were found independent of the size of f2, so here they hold independent of Ri. 
Using the Banach-Alaoglu Theorem we have the existence of a function 

« e L°°([0,T];L2(m")) nL2([o,r];i/i(M")) 

dtveL^{[Q,T]-{H^J{W^)) 

such that 

(8.12) v^^ ^vin L°°([0,T];L2(]R")) ^^ak* 

(8.13) v^' ^vin L^{[0, T]; i/^(M")) weakly 
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There exists an orthogonal basis {(pi} for i^([0, T]; (M")) where each function in the 
basis is smooth and has compact support in space. For Ri larger then the support 
of (j>, Theorem 18.41 with it's corollary show 



/ <v^\dt<j)> ds+ I < u"' ■ Vw^* ,(j)> ds+ [ <(j)- Vu^* , w^* > ds 
Jo Jo Jo 

T 

< \7v"' , V0 > ds =< vq, (j) > 

The limit to — )■ oo can be passed through the linear terms just as before. In 
the (compact) support of each basis function we have the strong convergence 
to pass the limit through the non-linear terms. A diagonal argument shows this 
convergence holds as i?j — > oo. □ 

In the above existence theorems, the pressure term can be found by either taking 
the divergence of the VCHE and solving the Poisson equation, or using a famous 
result of de Rham. See, for example, [28]. 



Theorem 8.7. The solutions to the VCHE constructed in Theorems \8.4\ and \8.6l 

with initial data in vq G , satisfy the bound 

(8.14) ||V*^t;||^+ f\\V''+'v\\l<C{n,a,,.,\\v\\^.) 

Jo 

for all M <K. 

Proof. We will do the calculations formally and note that these bounds can be ap- 
plied to the approximate solutions constructed in Theorem [821 this proof proceeds 
by induction. The inductive assumption is that the following bound holds for all 
TO < M. 

l|V'"f||2+ / \\V"'+\j\\ldt < c 
Jo 

The base case (to = 0) is true by Lemma [8.3[ we will now show that it holds for 
TO = AI. The bound [HU] with the inductive assumption implies 

f \\^"'+'u\\ldt<C 
Jo 

Multiply the VCHE ([TT]) by A^^w and integrate by parts to find 

(8.15) \±\\v^ty\\l + < Im + Jm 

hi = E O ^ ^"'^ ■ vv""'""^' v""'^ > 



m=0 
M 



jM = y { < v''v ■ vn, v''-"^v > 



m=0 



TO 



The two integrals on the RHS are estimated essentially the same way. The key 
difference is that in the first one we use the relation < u ■ Vw, v >= while in the 
second we can place an extra derivative on u. 
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With application of < w, Vu, v >— the first bound becomes 

/A A 



7n— 1 ^ ' 

Holder's inequality, the Sobolev inequality, and Cauchy's inequality show 

M 
M 

<C^\\V"^u\\l\\\7'''+'-^v\\l+'^\\\7' 

m—l 

Similarly for the second term 

M 



m=0 



M 



m=0 

Equation (|8.15p becomes 



M-m„,||2 



, M 

The Gronwall inequality with application of the inductive assumption finish the 
proof. □ 



Theorem 8.8. The solution to the VCHE constructed in Theorems \8.4\ and \8.6\ 

with initial data in vq G , satisfies the bounds 

(8.16) ||arV"t;||2 + ^*||arV™+iz;||^<C(n,a,^.,|l«|i^.) 

for all M + 2P <K. 

Proof. To prove this, we will bound the time derivatives of the solution in terms 
of the space derivatives, then use the previous theorem to establish regularity. We 
will do the calculations formally and note that these bounds can be applied to the 
approximate solutions constructed in Theorem 18.21 

Apply i9f V*^ to the solution of the VCHE, from this we have the inequality 



Using the Gagliardo-Nirenberg-Sobolev inequality and 18.11 we can bound the first 
term on the right hand side as 

P M 

< ^E E wdfv-^vrM-'y'^+'-^vwi 

p—O m— 
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Similarly for the second term, 



p—0 m— 
P M 

p—0 m—0 

Putting this together we can deduce 

^0 

This implies, for all M, P, such that M + 2P < K, 

\\dtv''v\\i<c\\vrH^. 

Appealing to Theorem 18 . 71 finishes the proof. □ 

The previous theorem demonstrates how the norms and ||w||^m+i can be 

bounded in terms of ||t'o||ff">- Since the PDE is parabolic we can expect regularity 
from interior estimates but the bounds will not depend explicitly on the initial 
conditions. 



Theorem 8.9. The solution to the VCHE constructed in Theorems \8.4\ and \8.6\ is 
unique. 



Proof. Let v and w be two solutions to the VCHE 11.11 with the same initial condi- 
tions. Let u and uj be the corresponding "filtered" velocities. The difference solves 
the PDE 

{v — w)t — v/\{v — w) + Vp + u ■ Vf — Lo ■ Vw + V ■ Vu'^ — w ■ 'S/uP" = 

with zero initial conditions. Multiplying this relation hy v — w and integrating by 
parts leaves 

^Wl + - "')ll2 =< {u-uj)- Vw, [v-w)> 

+ < {v — w) ■ 'S/{u — uj),v > 

+ < {v — w) ■ Voj, u — It; > 

Using Holder's inequality, the Gagliardo-Nirenberg-Sobolev inequality, Cauchy's 
inequality, and ()8.6|) . estimate the RHS 

< (u — uj) ■ Vw, {v — w) > < II M — aj||„|| Vu>||2||w — w|| _2ri_ 

<C\\v-w\\l\\Vw\\l+'^\\V{v~w)\\l 

< (v — w) ■V{u — uj),v > < \\v — w\\2\\V{u ~ aj)||„||t;|| 2,.^ 

<C\\v~w\\l\\Vv\\l+'^\\Viv--w)\\l 

< {v — w) ■ VUJ, V — W > < \\v — WII2II Vtilllnllw — It'll 



<C\\v^w\\l\\w\\l + -\\Viv-w)\\ 



2 



After using the bounds in Lemma 18.31 we have 



~\\v~w\\l+'^\\V{v-wm<C\\v-w\\l 
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By assumption ||wo — W0II2 = 0, so \\v — w\\2 — for all t S [0,T]. □ 
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